Condensate Fragmentation in a New Exactly Solvable Model for Confined Bosons 
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Based on Richardson's exact solution of the multi level pairing model we derive a new class of 
exactly solvable models for finite boson system. As an example we solve a particular hamiltonian 
which displays a transition to a fragmented condensate for repulsive pairing interaction. 
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Since the experimental realization of condensates of 
trapped bosonic atoms there exists a considerable grow- 
ing interest in refining the approximations to study fi- 
nite many body boson systems. Standard approaches 
used up to now are mean field based theories (MFT) like 
the Gross-Pitaevskii equation and the Bogoliubov theory 
which considers pair fluctuations on top of a Bose con- 
densate. In order to test such many body approximations 
and also to gain insight, exactly solvable models are ex- 
tremely useful. Two of such exactly solvable models were 
used in the context of boson traps, the one dimensional 
hard core boson model W and the harmonic interaction 
model 1^. Both model studies show the need to go be- 
yond MFT. We would like to present in this letter a new 
class of exactly solvable models based on a generaliza- 
tion of Richardson's exact solution of the pairing model 
[p|-p| . Of particular interest will be the model case with 
repulsive pairing interaction, where we find a fragmented 
condensate with two macroscopic eigenvalues of the one 
body density matrix corresponding to the lowest s and p 
wave states for bosons in a confining potential. The ocur- 
rence of condensate fragmentation in confined boson sys- 
tems is a much debated subject at present and so far no 
clear evidence for the existence of such a state has been 
demonstrated. Nozieres and Saint- James [|6| using the 
Hartree-Fock Bogolibov (HFB) approximation and later 
on Nozieres JT) using Hartree-Fock (HF) theory showed 
that fragmentation is energetically unfavoured with re- 
spect to a single condensate state for scalar boson sys- 
tems. Going beyond mean field Holzmann et.al. MM dis- 
cussed the possible existence a fragmented scalar boson 
condensate in two dimensions. The case of trapped spin 
1 boson systems has been treated in refs [^,Q where 
it was concluded that fragmentation might be realized 
in spin space but they still used mean field theory to 
describe the condensates in each spin component. 

Let us now come to the construction of our new ex- 
actly solvable model for interacing bosons. In rcf. |12] 
Richardson was able to obtain the exact eigenstates of a 
pairing hamiltonian for bosons 
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where ^A = J2a O-AaClAa , «A = J2a '^Aa'^Aa , (Aa) 

are the quantum numbers of a confining potential (in 



the particular example of a three dimensional spherical 
potential which we will discuss below A is the principal 
quantum number and a represents the orbital quantum 
numbers), and (Aa) is the time reversal of (Aa). The 
operator A]^ creates a zero angular momentum pair (sin- 
glet) in level A and the operator n\ counts the number 
of bosons in level A. 

In analogy with the solution for fcrmions Richardson 
assumed that the hamiltonian eigenstates can be written 
as a pair product wave function of the form 
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TV is the number of singlet 
boson pairs, Ei are pair energies to be determined by 
the eigenvalue equation and \ip) is a state of v unpaired 
bosons that fulfills Aa \(p) = and nA\ip) = t^Al^p). A 
given configuration in the Hilbert space is determined by 
the set VA of unpaired bosons in each level A. The total 
number of bosons is then A^ = 2M + J2a ^a- A crucial 
point is the observation that the off diagonal matrix ele- 
ments of the one body density matrix {'ip\ a\^aAia' |'0) ^^^ 
zero with |'0) given by yj^ because the action of a\^aA'a' 
on 1-0) either destroys a singlet pair increasing t/ by 2 or it 
acts on an unpaired particle chanching its configuration 
from (Aa) to (A'a'). In any case the state a\^aA'a' IV") 
with Aa 7^ A'a' will be orthogonal to [ip]. Therefore the 
occupation numbers {ip\ ha |V') Eire the eigenvalues of the 
one body density matrix. This situation is analogous to 
the Fermion case where this is known as the generalised 
seniority scheme |l3] . 

By acting with the hamiltonian on the trial wavefunc- 
tion, after a straightforward but long derivation, Richard- 
son obtained the following set of equation for the pair 
energies 
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where J7a is the degeneracy of level A .The energy eigen- 
values of (|l|) are given hy E = J^a^^i-'a + J2i=i ^i fo'^ 
each solution of the set of equations (H) as will be ex- 
plained below. 

The proof of integrability of the pairing model for 
fermion systems has been given by Cambiaggio et al. 



Q. They found a complete set of commuting opera- 
tors (constants of motion)in terms of which it is possible 
to express the pairing hamiltonian and the number op- 
erator as particular linear combinations. Unfortunately, 
since the authors were not aware of the Richardson's pre- 
vious works, they were not able to obtain the eigenval- 
ues of the new set of operators, nor the exact solution of 
the pairing hamiltonian. This connection has been estab- 
lished recently [n5| and the eigenvalues of the commuting 
operators found using Conformal Field Theory. 

The Richardson solution for fermions remained prac- 
tically unused till very recently when it was applied, in 
the context of ultrasmall superconducting grains |16|, to 
study the transition from the superconducting to the nor- 
mal state. We are not aware of any application for boson 
systems. 

In complete analogy to the fermion case, a set of global 
commuting operators can also be written in the boson 
case, based on the group algebra of the pseudo-spin gen- 
erators Kl = inA + {^K , K+ = IA{ = {K]^yoi 
5C/(1,1) . 

Ra=KI + 2.9 V i Ka ■ Ka, (4) 
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The scalar product in ( ^ refers to the SU{1, 1) group, 
Ka -Ka- = \ {K+Kl, + KIK+,) - K^Kl,. The set of 
operators Ra is complete, the operators commute among 
each other, and the pairing hamiltonian (|l|) can be writ- 
ten as the linear combination 2^^?7Ai?A. These con- 
ditions demonstrate that the pairing model is integrable 
for boson systems. 

We will now look for the eigenstates of the Ba opera- 
tors (Q) using the trial eigenstates (||) . The form of the 
boson pair amplitudes u\ in the boson pair operators 
Bi = J2a ""a^a is fixed by the solution of the one pair 
problem, namely u\ = 1/ (2riA — Ei). The eigenvalues 
^a IV") = -^a IV") can be worked out in an analogous way 
as in the original Richardson paper. The eigenvalues are 
given by 
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where Ea = ^a + 2i^A, while the pair energies Ei should 
be the roots of the Richardson equations (|^) with t^a re- 
placing £A- 

One can readily verify that for tja = £a in eqs. M) and 
(h), the linear combination 2 ^^ eaRa gives the Richard- 
son hamiltonian ffl) and 2 ^^ eaAa gives the correspond- 
ing eigenvalue. An important by-product of having found 
the eigenvalues (ra) is that the tja are in principle free pa- 
rameters not necessarily related to the single boson en- 
ergies. We will exploit this freedom to obtain solutions 
for generalized pairing hamiltonians. 



Looking back to the form of the hamiltonian (|i|) one 
sees that since the pair operators ^a are normalized to 
the square root of the degeneracy of the level A, the ef- 
fective pairing matrix elements are proportional to the 
square root of the product of the degeneracies of the two 
shells ^/TIaTIa'. In a spherical harmonic oscillator po- 
tential each degeneracy is in turn proportional to A'^"-'^ 
where d is the space dimension and A is the principal 
quantum number while the single boson energies are lin- 
ear in A, producing unphysical occupations of the high 
lying levels for attractive pairing or a compression of the 
bosons in the lowest level for repulsive pairing. We will 
obtain a more realistic hamiltonian making use of the 
freedom in the choice of the t^a. In order to cancel the 
undesired dependence of the effective pairing matrix ele- 
ments on the degeneracies we make the following defini- 
tion rjA = (ea)'^ in eqs. (Q) and (^. The new hamilto- 
nian, H = 2J2a^aRa, can be expanded using the defi- 
nition of the _Ra (§) in terms of the SU{1, 1) generators 
or equivalcntly in terms of the pair operators. The final 
form of the hamiltonian is (more details of its construc- 
tion will be given elsewhere) 
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with £ an uninteresting constant, ea = ea + 2Vaa ~ 
J2a' ^aa' ^A' being the single boson energies, and 

VAA'=9/2j:to^A4-'''- 

We can readily check that now the effective interac- 
tion terms in (O), contrary to (0), correctly scale with en- 
ergy. In addition to the pairing term in (0) also appears a 
particle-hole interaction of the monopole-monopole type 
which gives the hamiltonian a rather rich and quite gen- 
eral character. It has, however, the restriction that pair- 
ing and particle-hole interactions are linked to be of op- 
posite sign, a feature which may be realised only in par- 
ticular situations. On the other hand, we think, that the 
equality in magnitude of both interactions do not invali- 
date our general conclusions below. 

The energy eigenvalues of (m) can be obtained summing 
the eigenvalues (|]) as 
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It is worthwhile to emphasize at this point that the 
states (0) are common eigenstates of the operators R 
in (m and, consequently, of any linear combination of 
them like the pairing hamiltonian (Mj or our more gen- 
eral hamiltonian (0) provided that the pair energies are 
the solutions of equation (|). The solution of (ph was 
already discussed in ref. |l2 1 and we will give here only a 



brief summary. Assuming that we have L oscillator levels 
and 2J\f paired bosons, the M pair energies Ei are real 
roots of (||) in the interval —oo<Ei< 2rjL for an attrac- 
tive pairing interaction or 2?7o < Ei < oo for a repulsive 
pairing interaction. The groundstate has all pair energies 
in the restricted intervals —oo<Ei < 2r]Q for an attrac- 
tive pairing interaction or 2r]o < Ei < 2?7i for a repulsive 
pairing interaction. Any state in the Hilbert space cor- 
responds to a particular distribution of the Af pair ener- 
gies into the L intervals. States with broken pairs can be 
generated by replacing a boson pair B by two unpaired 
bosons which can occupy any of the L shells. For ex- 
ample, the first excited state with two unpaired bosons 
corresponds to solve (^ with Af — I pair energies in the 
first possible interval and ^0 = 1^1 = 1. 

Having in mind that the eigenstates of the hamiltonian 
(^ are the same as those of the pairing hamiltonian (|^), 
the occupation numbers for a given state can be calcu- 
lated as the derivatives of the pairing hamiltonian Hp 
(0) with respect to the the single boson energies e as has 
been done by Richardson in Q 
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Details of the derivation of the occupation numbers 
and the final set of equations can be found in ref. |12[ . 
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FIG. 1. Occupation numbers for 1000 bosons in 50 har- 
monic oscillator shells, the interaction strength is g = —2.0, 
giving a depletion factor of 0.608. 

We have performed a series of calculations for a system 
of 1000 bosons {N = 500) trapped in a three dimensional 
spherical harmonic oscillator with a cutoff at 101/2 fiu) 
(L = 50). For attractive pairing all quantities evolve 
smoothly with increasing attraction. In fig. 1 we show 



the occupation numbers for g — —2.0 corresponding to 
a depletion factor of 0.608. They display a reasonable 
pattern filling first the lowest levels. This is not the case 
for the pairing hamiltonian (nf) with t^a = £a which for 
d > 1 populates first the high lying levels. The compar- 
ison between the exact results and approximations like 
Hartree-Fock-Bogoliubov or their number conserving ex- 
tensions |jl^,|l^ will be the topic of a future work [^ . 

Changing now to repulsive pairing we have found an 
unexpected feature. For some critical value of g the nor- 
mal groundstate boson condensate suddenly turns into a 
state in which the bosons are condensed into the A = 
and the A = 1 states while the occupation of the other 
levels is negligible. This state was already envisaged by 
Richardson himself studying an approximate solution for 
eq. (S) in the thermodynamic limit |13] . Taking into ac- 
count that the one body density matrix is diagonal in the 
basis (Aa) (see above) the occupation numbers in (^) are 
their diagonal matrix elements or, equivalently, the den- 
sity matrix eigenvalues. Therefore, this new state with a 
macroscopic occupation of the two lowest harmonic os- 
cillator shells constitutes a truly fragmented condensate 
state. It is commonly accepted since the work of Nozieres 
and Saint James g] that for homogeneous systems frag- 
mentation cannot occur in systems of scalar bosons with 
repulsive interactions. This might perhaps be the first 
example of fragmentation in a confined boson systems. 

In order to understand the nature of this new 
quantum phase we consider a coherent state |(/)) = 

exp yj2AfT^ |0) where F is the most general time rever- 
sal invariant coherent boson that breaks rotational and 
reflection symmetry. 
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where the three independent variables are defined in 
the intervals /3>0, 0<7<7r and 0< Lp < 2tt. 

Since, as mentioned before, the groundstate is a com- 
mon eigenstate to our hamiltonian (m and to the pairing 
hamiltonian (m) with t^a = {hcoA + 3/2) , for simplicity 
we minimize the energy of the latter hamiltonian which, 
indeed, will turn out to be a very accurate procedure. 
Apart from a constant energy term, the groundstate en- 
ergy and the occupation numbers are given by 
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where we have defined the adimensional parameter x = 
^^. The energy ( pl3| ) is independent of the variational 
parameters 7 and (^. Minimization of (nOh with respect 

to (3 gives the solutions [3 — for x <1 and (3 — y fry 
for a; > 1. The critical interaction strength correspond- 
ing to the phase transition at a; = 1 is ^c = f^- On 
the other hand, it is easy to check that the Bogoliubov 
approximation has a break down at this critical value of 
g showing the instability of the singlet boson condensate 
against pair fluctuations. Beyond the critical point the 
occupation numbers for the first two levels, obtained by 
inserting the minimum (3 value in ( |ll| ) , are ng 
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FIG. 2. Occupation numbers no and n\ for 1000 bosons in 
50 harmonic oscillator shells as a function of the adimensional 
parameter x. 

In figure 2 we have plotted the occupation numbers of 
the first two levels for a system of 1000 bosons interact- 
ing in 50 levels of a three dimensional harmonic oscillator 
as a function of the parameter x . The exact results in 
solid lines clearly display an abrupt transition for x — 1 
as predicted by the meanfield description. In fact the ap- 
proximate results are indistinguishable from exact ones. 
Moreover the total occupation of the rest of the levels is 
always lower than 10~^. We have added two insets to the 
figure to show a close up of the transition region for the 
occupations uq and ni. The small differences between 
the exact and the approximate results suggest that this 
is a 1/J\f effect, and that we are seeing the precursor of 
a true quantum phase transition. 

The new phase is a rather peculiar fragmented state 
characterized by a macroscopic occupation of only the 
two lower levels. It can be approximated by a single 
condensate state at the price of breaking of the refiection 
symmetry generating a permanent dipole deformed state. 
As we have seen before the parameter x is proportional to 
the number of bosons and inversely proportional to the 



oscillator frecuency. Assuming that this phase transition 
might be realized in realistic systems, it can be controlled 
by varying the number of bosons or the characteristics of 
the confining potential. 

In conclusion, we have developed a new class of exactly 
solvable models for confined boson systems. These mod- 
els are exactly solvable in any dimension (refiected in the 
degeneracies Vl/y) and with any kind of confining poten- 
tial (reflected in the single boson energies ea). We have 
made numerical applications showing the great potential 
utility of the model for serving as a testing tool for many 
body approximations. Moreover we have found a quan- 
tum phase transition to a fragmented state for repulsive 
pairing interactions. To the best of our knowledge, this 
is the first example of fragmentation in a confined scalar 
boson system. As such, it may stimulate further investi- 
gations to see whether this new phase can arise in more 
realistic situations. 
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